UNIMODULARITY AND PRESERVATION OF VOLUMES IN NONHOLONOMIC 

MECHANICS 
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Abstract. The equations of motion of a mechanical system subjected to nonholonomic linear constraints can be 
formulated in terms of a linear almost Poisson structure in a vector bundle. We study the existence of invariant 
measures for the system in terms of the unimodularity of this structure. In the presence of symmetries, our approach 
allows us to give necessary and sufficient conditions for the existence of an invariant volume, that unify and improve 
"i results existing in the literature. We present an algorithm to study the existence of a smooth invariant volume for 

g nonholonomic mechanical systems with symmetry and we apply . to several concrete mechanical examples. 
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1. Introduction 

The existence of an invariant measure for a system of differential equations is a very important 
property. From the point of view of dynamical systems, it is a key ingredient for the application 
of ergodic theory. It is also a crucial hypothesis in Jacobi's theorem of the last multiplier that 
establishes integrability of the system via quadratures, see e.g. [2]. Moreover, the existence of a 
smooth invariant measure imposes certain restrictions on the qualitative nature of the fixed points of 
the system; namely, it prohibits the existence of asymptotic equilibria and limit cycles. The methods 
developed in this paper address the general question of the existence of an invariant measure for 
nonholonomic mechanical systems. 

Suppose that a given system of differential equations is described by a vector field X on an 
orientable phase-space manifold M that is equipped with a volume form fi. If / is a positive function 
on M, then X preserves the measure ffi if and only if / satisfies Liouville's equation 

/ + /div^(X) = 0, (1.1) 

where / = X{f) is the derivative of / along the flow of X. In local coordinates, the above equation 
is a linear partial differential equation for /. We remark that the existence of an invariant measure is 
a global property of the system. By the Flow-Box Theorem, Liouville equation can always be solved 
locally (away from equilibrium points). 

Liouville equation for nonholonomic systems with symmetry 

Let us summarize the local version of our results. In the case of a nonholonomic mechanical system 
with linear constraints, the phase-space manifold M is a vector bundle D* over a configuration 
manifold Q. The fiber variables p are physically interpreted as momenta. 

For simplicity, we restrict our attention to systems without potential forces. We are especially 
interested in systems that possess a symmetry group G that acts on Q leaving the constraints and 
the kinetic energy function invariant. In section |4] we carry out the details of the reduction of the 
equations of motion. It turns out that the reduced phase-space D* /G is a vector bundle over the 
shape space Q := Q/G, and the fiber variables can be split into (payPa) where, roughly speaking, the 
Pa are compatible with the symmetries. 

The Hamiltonian H of the system is derived from the kinetic energy and is a positive definite 
quadratic form, with our choice of {pa,Pa), is of block type H = ^{9°''^{(l'')PaPi3 + 9°'''{(j'')PaPb), where 
g'' are local coordinates on Q. The reduced equations of motion that obey the Lagrange-D'Alembert 
principle have the form 

dt "9p„' 

dPa__pdH J dH 
dt - ^'^dq^ ^'^'P'dpj' ^ ^ ^ 
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where the upper-case indices I, J, . . . run over the joint range of a, 6, . . . and a, (3, ... , while Y^ and 
Cfj are smooth functions of g^ Moreover, the functions Cfj are skew-symmetric with respect to the 
lower indices, i.e. Cfj = —Cjj. 

The Lagrangian version of the previous equations was considered in [26] (see also [16]) for the 
reduction of a nonholonomic mechanical system. In the particular case when the indices t run over 
the same range that the indices a, that is, if dimension assumption holds (see [1]) then the previous 
equations were considered in [2] 



As we shall see (consequence of Theorems 3.6 and 3.8), if there is a preserved measure $ for the 



system (1.2) then it may be chosen of basic type, i.e. 

$ = e^^"^') dq' A dpa A dpa, 
for a certain real smooth function a on Q that needs to be determined. If we enforce the Liouville 



equation (1.1) (for / = e°", /i = $ and X the vector field defined by equations (1.2[)), we get 



SY^SH ,^^^aH_ 




d(f dpa ° dcf dpa 

where we have used the skew-symmetry of the coefficients Cfj and the equality of mixed partial 
derivatives of H. Rearranging terms we get 

dq^ + ^'^dq^ + ^^') dp^ + ^^'dpa ~ 

The form of the Hamiltonian H implies that the latter identity can only hold for all values oi Pa,Pa 
if the following identities hold separately: 

BY'- -^ Bd 

-^ + ya^^^+Clr = () for all a, and Cf . = for all a. 

oq'- oq'- 

This paper is devoted at giving an intrinsic version of the above conditions and applying them to 
several concrete mechanical examples. Notice that only the first condition involves the unknown a 
and that none of them involves the Hamiltonian function H. All of the relevant information for the 
existence of an invariant measure is contained in the coefficients Y^ and Cfj. We shall see that these 
coefficients determine a linear almost Poisson structure on D* /G. The formulation of nonholonomic 
mechanics in terms of a linear almost Poisson structure is well-known and has its origins in [33j (see 
also [9l UHl 1211 EH]). In this paper we show that the above conditions for a preserved measure can 



be formulated intrinsically (Theorem 4.3) in terms of a geometric notion: the unimodularity of the 
underlying linear almost Poisson structure. 

Unimodularity and measure preservation in mechanics 

The well known theorem of Liouville states that a Hamiltonian system on a symplectic manifold 
preserves the symplectic volume, see e.g. [1]. The situation is not so simple for a Hamiltonian system 
on a Poisson manifold. For instance, if the Poisson manifold is the dual space q* of a Lie algebra q 
equipped with the Lie-Poisson structure, then Kozlov ^25j showed that the flow of a Hamiltonian of 
kinetic type on q* preserves a smooth measure if and only if the Lie algebra q is unimodular. 

More generally, a sufficient condition for the existence of a smooth measure for a Hamiltonian 
system on an abstract Poisson manifold M is that the modular class of M vanishes, see e.g. [39] . 
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The modular class is the Poisson cohomology class of the modular vector field of M with respect 
to a volume form. This vector field plays an important role in the classification of certain Poisson 
structures (see [121 113 1211 )• On the other hand, the modular class of M is the obstruction to 
the existence of a duality between the canonical homology and the Poisson cohomology of M (see 

The modular class of a Lie algebroid was introduced in p^J and generalizes the modular character 
of a Lie algebra. In Marrero [29] the results of Kozlov [25] are generalized to consider the preservation 
of volumes for mechanical Hamiltonian systems on vector bundles equipped with a linear Poisson 
structure, which is the dual object of a Lie algebroid. A Hamiltonian is mechanical if it can be 
expressed as the sum of the kinetic and potential energies. Just as in the case studied by Kozlov, 
the unimodularity of the Lie algebroid is intimately related with necessary and sufficient conditions 
for the existence of an invariant measure. 

Contributions and outline of the paper 

In section [2] we introduce the notion of modular vector field with respect to a volume form and 
unimodularity of a general almost Poisson structure. An almost Poisson structure on a manifold M 
is unimodular if there exists a volume form on M which is preserved by all the Hamiltonian vector 
fields. In section |3] we specialize to the case when the almost Poisson structure is linear in a vector 
bundle and we analyze the relationship between preserved measures for Hamiltonian vector fields 



and unimodularity. The main theoretical result of the paper is contained in Theorem |3.8[ that shows 
that a kinetic energy Hamiltonian vector field on a vector bundle possesses an invariant measure if 
an only if the underlying almost Poisson structure is unimodular. Recently, Grabowski [T3] proved, 
independently and using a different approach, a similar result. He introduced the notion of the 
modular section of a skew-symmetric algebroid and from this theory he obtained the result. 

Section |4] focuses on preservation of volumes for nonholonomic systems. After carrying out the 
reduction by a symmetry group from the almost Poisson perspective, we state and prove Theorem 



4.3 and Corollary 4.4 that allow us to recover several known results in the literature [SI [201 ESI 112] in 
a unified framework. Moreover, we are able to express the related necessary and sufficient conditions 
for the existence of an invariant measure both in global and local form, and to extend the results to 
the case where the dimension assumption (see [1]) does not hold. 

In section [5] we give an algorithm to study the existence of an invariant measure for the reduction 
of a nonholonomic system with symmetry, that is applied in section [6] to several concrete mechanical 
examples. This allows us to obtain the following original results: 

We prove non-existence of an invariant measure for the reduced equations of motion of the following 
problems: 

(i) As considered in [371 138]. the motion of a rigid body with a planar section, that is not planar 



or axially symmetric, and rolls without slipping over a fixed sphere (Theorem 6.1). 
(ii) The motion of an inhomogeneous sphere whose center of mass does not coincide with its 
geometric center that rolls without slipping on the plane (Chaplygin's top) and is not axially 



symmetric (Theorem 6.2). 



(iii) The rolling without slipping of an inhomogeneous sphere, whose center of mass coincides with 



its geometric center, over a circular cylinder (Theorem 6.3). 



(iv) We also show that the measure fails to exist in the limit when the radius of the cylinder in 



(iii) goes to zero and the given ball is axially symmetric but not homogeneous (Theorem 6.5). 
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This example utilizes our extended theoretical results since the dimension assumption does 
not hold. 

On the other hand, we consider the motion of an axially symmetric Chaplygin top subjected to 
an additional rubber constraint and we give an explicit formula for the preserved measure (formula 



(6.13)). 



2. Almost Poisson structures and unimodularity 

We begin by introducing the notion of the modular vector field of an almost Poisson structure with 
respect to a volume form. This is the first step towards the definition of the concept of unimodularity. 

Recall that an almost Poisson structure on a manifold M is an M-bilinear bracket of functions, 

{-, ■} : C^{M) X C^{M) -^ C°°(M), 

that satisfies 

(i) skew-symmetry: {f,g} = -{fi-, /}, 
(ii) Leibniz rule: {fh,g} = f{h,g} + h{f,g}, 

for all f,g,h G C°^{M). We speak of an almost Poisson bracket since we do not require the Jacobi 
identity 

{{f,9},h} + {{h,f},g} + {{g,h},f} = 

to hold. If the above identity holds for all functions, then {■, ■} is a usual Poisson structure on M 
|35j . Our interest in almost Poisson structures comes from nonholonomic mechanics ^33] (see also 

As a consequence of Leibniz' rule, for any function / e C°°{M) we can associate a vector field Xf 
on M, the Hamiltonian vector fielcn of f, whose action on g E C°°{M) is defined by Xf(g) := {g, /}. 
Note that by skew-symmetry we have Xf{g) = —Xg{f). 

Let (x*) be a local system of coordinates on M, and denote by vrjj(x) := {x'^,x^}. The local 
expression for Xf is 

and the integral curves of Xf satisfy Hamilton's equations 

CLX . , OJ . , 

- = vr.,(x) — (X). 

Now suppose that M is an orientable manifold and that $ is a volume form on M . The following 
proposition is key to the developments in this paper. 

Proposition 2.1. The map M$ : C°°{M) -^ C°^{M) given by 

U^{f) = div$(A^) 

defines a vector field on M, where div$(Xj) denotes the divergence of the vector field X f with respect 
to the volume form $. 



Note that Xf is not a Hamiltonian vector field in the usual sense since we work with an almost Poisson bracket. It would be more 
accurate to talk about an almost Hamiltonian vector field. We eliminate the "almost" in our terminology for brevity. 
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Proof. We need to show that M$ is a derivation, i.e. that 

MM9) = mM+9y^M)- (2.1) 

We shall see that this is a consequence of the skew-symmetry of the bracket. The reader can check 
that the above identity holds with a simple calculation in local coordinates. Here we present an 
intrinsic proof. Denote by £xC( the Lie derivative of the form a along the vector field X. Using the 
properties of the Lie derivative we have: 

M$(/^)<l> = £x,,<^ = £fx,<^ + £,x,^ 

= {fMM + g^M)) <l> + dfAix,<^ + dgA ix,$. ^ ■ ' 

Using that (i/ A $ = 0, one obtains 

= ixMf A <^) = X,{f)<^ -dfA ix,$, 
which implies 

dfAix,<^=Xg{f)<l>. 



Similarly, dg A ix,^ = Xf{g)^. Substituting the latter expressions into (2.2) and using the skew- 



symmetry of the bracket yields (2.1). D 



If dimM = m and the system of coordinates (x*) is such that $ = dx^ A ■ ■ ■ A dx"^, then the local 
expression for M$ is 

Definition 2.2. The vector field M$ defined above is called the modular vector field of the almost 
Poisson structure on M with respect to the volume form $ . 

Remark 2.3. The notion of modular vector field for a Poisson manifold has been extensively used 
J12t \rj\ [27] . Moreover, Weinstein shows in [32] that it defines a Poisson cohomology class of order 1 
which is independent of the choice of volume form. The Poisson manifold is said to be unimodular if 
this cohomology class is zero, or, equivalently, if the modular vector field is Hamiltonian. o 

It is natural to ask how the modular vector field is affected when we consider a different volume 
form $' = 6*^$ on M . Here a is any smooth function on M . A straightforward calculation gives 

M$/ = Me^$ = M$ - X„. (2.4) 



The above calculation, together with the comments in Remark 2^, suggest that we make the 
following definition: 

Definition 2.4. The orientable almost Poisson manifold M is said to be unimodular if there exists 
a volume form ^ on M and a real function o G C°°(M) such that M$ = X^j. 

The following Theorem is an immediate consequence of the definition. 

Theorem 2.5. // the orientable almost Poisson manifold M is unimodular there exists a volume 
form on M that is invariant by the flow of all Hamiltonian vector fields. 

Proof If M$ = X^, then diVe<.$(X/) = for any function / G C°^{M). D 
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3. Linear almost Poisson structures, kinetic energy hamiltonians and 

unimodularity 



Theorem 2.5 shows that unimodularity is a sufficient condition for the existence of an invariant 
measure for Hamiltonian vector fields. In this section we investigate the degree to which this condition 
is also necessary in the special case when we deal with linear almost Poisson structures on vector 



bundles. We will prove in Theorem 3^ that if the Hamiltonian vector field associated to a kinetic 
energy Hamiltonian preserves a smooth measure then the corresponding almost Poisson structure is 
unimodular. This is the main theoretical result of the the paper and will be applied in section [6] 
to analyze the existence of an invariant measure of several nonholonomic mechanical systems with 
symmetry. 

3.1. Linear almost Poisson structures and unimodularity. Let t : E ^ Q he a. real vector 
bundle of rank n on a manifold Q of dimension m. We shall denote by Eg the fiber over the point 
q E Q and by r(r) the space of sections oi r : E ^ Q. 

A real function if : E ^ M. is said to be linear if its fiberwise restriction (p\^ : i^g — )■ M is linear 
for all q ^ Q. On the other hand, a real function 5" : i? — )■ M is said to be basic if it can be written 
as o" = (J o r for a certain a G C°°{Q). 

It is well known that there exists a one-to-one correspondence between the space of sections T{t*) 
of the dual bundle r* : E* ^ Q and the space of linear functions on E. Such a correspondence is 
the following: to the section X G r(T*), we associate the linear function X^ : E' — )■ M defined by 

Xe{a) = X(r(a))(a), for a e E. 

For every point q E Q we can choose local coordinates (g*) on a chart U and a local basis {e/} of 
r(r*) on U . We denote by p/ the local linear functions: 

Pi = {ei)i. 

Then {q^,pi) is a system of local bundle coordinates on (r*)^^(t/) C E, where g* and pi are basic 
and linear functions, respectively. 

Definition 3.1. ^26j An almost Poisson structure on a real vector bundle r : E ^ Q is said to be 
linear if the almost Poisson bracket of linear functions is linear. 

The following proposition explains how basic and linear functions behave under the action of linear 
almost Poisson brackets. Its proof is given in 



Proposition 3.2. [26] Let {-, ■} be a linear almost Poisson bracket on E. 

(i) If X is a section of t* : E* -^ Q and a is a smooth function on Q, then the bracket {cror, X^} 
is a basic function on E (the bracket of a linear function and a basic function is a basic 
function). 

(ii) // cr, 0"' are smooth functions on Q then {a o r, a' o r} = (the bracket of basic functions is 
zero). 



The content of the above proposition can be stated in terms of the local bundle coordinates {q\ pi) 



as 



{pi,Pj} = -CfjPK, {q\pi} = p], {q\q'} = 0, (3.1) 
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for certain smooth functions Cfj, p\ on Q. These functions are called the local structure functions 
of the linear almost Poisson structure. In view of this, we obtain the following expression for the 
Hamiltonian vector field X^ of a smooth function H : E ^ M.: 

and the corresponding Hamilton's equations 

dq' jdH dpi f ^dH , ^^ dH 



It is hard to underestimate the importance of the above universal form of equations in mechanics. 
The equations of motion of a great number of mechanical systems can be formulated as above. 

Let a G C°°{Q). For further reference we note that the Hamiltonian vector field of the basic 
function cr o r : i? —>■ M is the vertical vector field 

Basic-unimodularity 

Suppose that the manifold Q and the vector bundle t : E ^ Q are orientable (then the same is 
true about the dual bundle t* : E* ^ M.). Denote by «/ a volume form on Q, and by fi G r(A"r*) 
a volume form on the fibers of E* (that is, Vt is an everywhere non-vanishing section of the vector 
bundle A"r* : hJ^E* — )■ Q, where n is the rank of E). As indicated in the Appendix, we can consider 
a basic volume form $ = z/ A $7 on £" (actually any basic volume form on E can be obtained through 
this construction). 

Local expressions for the modular vector field M<j, are obtained as follows. Recall that we had 
selected a chart U with local coordinates (g*), and that {e/} was a local basis of r(r*) on U . Suppose 
moreover that locally we have 

V = dq^ ^■ ■ ■ ^ dq"^, fi = Ci A ■ ■ ■ A e„. 

Then (see Appendix), 

$ = z/Af2 = rfg^A---A dq"" A dpi A ■ ■ ■ A dp^- 



Under these assumptions, the local expression (2.3) for the modular vector field M<j, becomes 



^""-if + ^")|;- P-^) 



As a consequence we have 



Proposition 3.3. The modular vector field of a linear almost Poisson structure on a vector bundle 
T : E -^ Q with respect to a basic volume form is r-vertical. 

The following definition specializes our interest in the relationship between basic volumes and 
unimodularity. 

Definition 3.4. Consider a vector bundle r : E ^ Q equipped with a linear almost Poisson structure. 
The structure is said to be basic-unimodular if there exists a basic volume form $ on E, and a basic 
function (TOT eC^iE) (with a G C^iQ)) such that M$ = X„ 



»-o"or • 
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Remark 3.5. Suppose that a linear almost Poisson structure on a vector bundle t : E -^ Q is 



basic-unimodular and that $' is an arbitrary basic volume form on E. Then, from (2.4), it follows 
that there exists a' G C°°{Q) such that M$/ = X„ 



^cr'OT- 



The following theorem shows that unimodularity, basic-unimodularity, and the fact that all Hamil- 
tonian vector fields preserve a basic volume are equivalent notions for linear almost Poisson structures 
on vector bundles. 

Theorem 3.6. Let t : E ^ Q be a vector bundle equipped with a linear almost Poisson structure. 
Suppose that both Q and E are orientable. Then the following statements are equivalent: 

(i) The almost Poisson structure is unimodular. 
(ii) The almost Poisson structure is basic-unimodular. 

(iii) There exists a basic volume form that is invariant under the action of all the Hamiltonian 
vector fields. 

Proof. The implications {Hi) =^ (ii) and (ii) =^ (i) are trivial. Moreover, (ii) =^ {Hi) is a 



direct consequence of Theorem 2.5 Therefore, it only remains to prove the implication (i) =^ (ii) 



Suppose that the linear almost Poisson structure on E is unimodular. Again, by Theorem 2.5[ 



there exists a volume form ^ on E that is preserved by all Hamiltonian vector fields. Such form can 
be written as 

$ = e^u AQ, 
where i/ is a volume form on Q, Q ^ r(A"r*) is a volume form on the fibers oi t* : E* ^ Q, and a 
is a smooth function on E. We shall prove that a can be chosen to be basic. First note that for any 
Hamiltonian function / G C'^{E) we have 

= (div$X;)<l> = i?x,$ = {Xf{d) + M,Af7(/)) $. 

Hence, X^{f) = "Mu^nif) for any function / G C°°{E) which implies 

Xa = M,^n. (3.6) 



It follows from Proposition |3.3| that the vector field X^ is r-vertical. Now let : Q — ?■ £" be the 
zero section of the vector bundle r : E ^ Q. Consider the function a := 5" o G C°°{Q). Using 
the verticality of X^ together with the local expression for the Hamiltonian vector fields (|3.2) one 



obtains X^ o = X^ot ° 0. Therefore, evaluation of (3.6) along the zero section yields 



Since both Xo-or and M^/^^ are vertical lifts of vector fields in Q, the above identity implies X^, 
Mi^AH- In other words, the almost Poisson structure is basic-unimodular. 



D 



The above theorem, together with the local expressions (3.4) and (3.5), give us the following local 
characterization of unimodularity. 

Corollary 3.7. A linear almost Poisson structure on an orientable vector bundle r : E ^ Q is 
unimodular if and only if there exists a smooth function a : Q — > M such that for any system of 
bundle coordinates {q\pi) we have 

''W^t^^^"^^' ^''''" / = l,...,n. (3.7) 
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3.2. Kinetic energy Hamiltonians and unimodularity. Now suppose that the vector bundle 
r : -E — )■ Q is also endowed with a bundle metric S '■ E x q E -^ M.. For each g G Q, we have a scalar 
product S(g) : Eg x Eg -^ M.. In applications to nonholonomic mechanics this metric comes from the 
kinetic energy of the problem. 

The bundle metric S defines a Hamiltonian function H on E. Namely, the kinetic energy associated 
to 9 that is given by H{a) = ^9{q){a, a) for a G Eg. In terms of our local bundle coordinates {q\pi) 
we have 

H{q\pi) = Is'"" {q')pjPK, (3.8) 

where the g-dependent coefficients S'^^(g*) are obtained from the local basis {e^} of r(r) as 5'^^{q) : = 
S(g)(e'^(g), e^{q)). Note that the condition that S is a bundle metric implies that the matrix S'^^{q) 
is symmetric and positive definite for all q & Q. In particular, it is non-degenerate. 

We are now ready to state and prove the main theoretical result of the paper that relates the 
existence of an invariant measure for the Hamiltonian vector field of H defined as above, and the 
unimodularity of the underlying linear almost Poisson structure. 

Theorem 3.8. Let r : E ^ Q he a vector bundle equipped with a linear almost Poisson structure. 
Suppose that both the vector bundle and Q are orientable. Let S be a bundle metric on E and let 
H : E -^ M. be the kinetic energy Hamiltonian associated to S- Then the following conditions are 
equivalent: 

(i) The almost Poisson structure is unimodular. 

(ii) There exists a volume form ^ on E that is invariant under the flow of the Hamiltonian vector 
field Xfj . 



Proof. That (i) =^ (ii) follows from Theorem 2.5 To show the converse fix a volume form z/ 



in Q and a volume form fi G r(A"'r*) on the fibers of E*. The volume form $ can be written as 
$ = e'^u A i7 for a certain function a G C°°{E). Choose local coordinates (g*) on Q and a local basis 
{e/} of r(r*) satisfying 

z/ = (ig^ A ■ ■ ■ A rfg™, fi = d A • • ■ A e„. 

Denote by {q\pi) the corresponding system of fiber coordinates on E. We then get 

$ = e^'dq^ A ■ ■ ■ A dg™ A dpi A ■ ■ ■ A dp^. 



Considering the above expression for $, together with the local expression (3.2) for the Hamiltonian 
vector field X^, then the condition div$(X//) = yields 



dpiX'dqi dq^ 'V dpiX'dq^ '''^'"dp 

Differentiating this expression with respect to pl gives 
d^H ( ^ do dp\ ^j\ dE ( , d^a 



dpidpi \ dq^ dq^ J dpi \ dp^dq 

" dpLdpi Y'dq^ ^^"^""Opj) + dp J Y' dpLdt^^'' Opj^^''^"" dpj^dp. 



(3.9) 
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We now use our hypothesis that the Hainiltonian H is of the form (3.8). Omitting the dependence 
on g* we have 

d^H ^Lj dH _ J J dH 1 (93^-^ 



S^% -7^ = 9''PJ, -7^ = ^^^PIPJ. 



dpidpi dp I dq^ 2 dq 



Next we evaluate (3.9) along the zero section : Q ^ E and use the above relations to obtain 



^Li I ^ I 9a ^\ , dp} J 



Consider the smooth function o" on Q defined as a = a o 0. Then — -^ o = 7-^. Therefore, using 

oq'^ aq^ 

the non-degeneracy of S^^, we obtain 

p}|^ + 1^ + C/j = for all / = !,..., n 



(9g* dq^ 



, .1,, 



which is the local characterization for unimodularity given in (3.7). D 



4. Invariant measures for nonholonomic mechanical systems with symmetry 

In this section we will layout the almost Poisson reduction of nonholonomic mechanical systems 
with symmetry. We are particularly careful in the development of the associated geometry, to 



formulate in Theorem 4.3 the necessary and sufficient conditions for the unimodularity of the re- 
duced almost Poisson structure. This is the key step for our study of preservation of volumes for 
nonholonomic mechanical systems. Our result allows us to recover the existing results in the area 
[HI EQI EHl II2] in a unified framework, and to further extend them to the case where the dimension 
assumption does not hold. 

4.1. Reduction of nonholonomic mechanical systems with symmetry. We begin by intro- 
ducing our geometric definition of a symmetric nonholonomic system. 

Definition 4.1. Let p : Q ^ Q := Q/G be a principal G-bundle. A symmetric nonholonomic 
system is a couple (S,-D), where S is a G-invariant Riemannian metric on Q and D is a G-invariant 
constraint distribution on Q. 

Notice that the theory for a general nonholonomic system without symmetries is recovered by 
supposing that the symmetry group is trivial G = {e}. We will describe the dynamics of the system 
for a general symmetry group G in terms of linear almost Poisson structures. As mentioned before, 
this formulation goes back to [33] (see also [SI UHl 1211 EH])- 

The generalized nonholonomic connection 

For our purposes it will be convenient to have a clearer understanding of the structure of the 
constraint space D and its relation with the infinitesimal symmetries defined by the group action. 
Denote by Vp the vertical subbundle of p. That is, 

Vpiq):=T,OrhGiq), 

where Oihciq) denotes the orbit of G through q E Q. We will assume that the intersection 

V^p:=VpnD 
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has constant rank along Q. In addition to this hypothesis, a number of important works (e.g. |1], 
(i2] ) assume that the dimension assumption holds, namely that 

Vp + D = TQ. (4.1) 

In our approach we will not assume that the latter condition holds in general. 

Let CK := {V^p)-^ fl D be the orthogonal complement of V^p in D. Under our assumptions, CK is 
a G-invariant subbbundle of D that allows us to decompose 

D = V^p® J{. (4.2) 

We define a principal connection, the generalized nonholonomic connection, whose horizontal dis- 
tribution is given by 

S):='K®{yp®'K)^. 

In fact, since TQ = (Vp © :K) © (Vp © :K)^, it follows that 

TQ = Vp®Sj, 

and, using the G-invariance of the bundle metric S, one also checks that the horizontal distribution 
Sj is G-invariant. 



If the dimension assumption (4.1 ) holds, then the horizontal distribution S) = "K, and the principal 



connection defined above coincides with the the nonholonomic connection introduced in |1] . 



By G-invariance of the vector bundles Vp and !K, the decomposition (4.2) drops to the quotient 

D:=D/G = V^p/G®^/G. (4.3) 

The space V^p/G is a vector subbundle (over Q) of D that we denote by r^Op/Q : Vp/G — )■ Q. 
The space of sections r(r-y-Dp/G') is naturally identified with the G-invariant vector fields Z on Q that 
lie on the intersection V^p = D n Vp. As a consequence, for any Z G r(rvOp/G), the projection Z 
onto Q via p vanishes. This fact will be used in the proof of Theorem 



4.3 



below. 



Now consider the distribution "K on Q whose characteristic space at the point q & Q is 

3<{q) = {Tqp){'K{q)), with g G Q and p{q) = q. 
The above definition of 'K makes sense since "K is G-invariant. Notice that if the dimension assump- 



tion (|4^ holds, then 'K{q) = TgQ. 

On the other hand, the linear map (Tgp)|j{(g) : 3i(g) — ?■ H{q) is an isomorphism of real vector 
spaces. In fact, the horizontal lift ^ : 0-C{q) — )■ !K(g) associated to the generalized nonholonomic 
connection is just the inverse morphism of (Tgp)|iK(g). Thus, the horizontal lift of a vector field Y on 
Q that lies on CK is a G-invariant vector field Y^ on Q which belongs to the distribution [K. That 
is, Y^ G T^Tji/c), where T'k/g '■ ^/G — )■ Q is the vector bundle projection. 

Therefore, the vector bundles "K/G -^ Q = Q/G and "K ^ Q = Q/G are isomorphic and the 
tangent map to p, Tp, induces a vector bundle isomorphism between them. 
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Construction of the reduced almost Poisson bracket 

It is well-known that the space of orbits T*Q/G of the cotangent lifted action of G to T*Q is a 
vector bundle over Q that admits a linear Poisson structure {•, ■}t*q/Gi see e.g. [31]. We shall denote 
the natural bundle projection by f* : T*Q /G — )• Q. 

Note that the space of sections of the dual bundle f : TQ/G -^ Q may be identified with the space 
of G-invariant vector fields on Q. Let X and Y be any such vector fields and f,k & G°°{Q). Then 
the bracket {■, -yr'Q/G is characterized by the relations 

{Xe,Y,}T*Q/G = - {[X,Y]), , {fof*,Xe}T'Q/G = X{f)of*, {f of* ,kof*}T*Q;G = 0, (4.4) 

where X denotes the projection of the G-invariant vector field X onto Q via p. In the case where 
the Lie group G is trivial, the above formulas characterize the canonical Poisson structure in T*Q. 

Now denote by t* : T*Q — t- D* and J"* : D* — t- T*Q the dual morphisms of the canonical inclusion 
L : D "-^ TQ, and the S-orthogonal projector 7 : TQ — > D. Since D and S are G-invariant, 
the cotangent lifted action to T*Q induces a G- action on D* and the morphisms l* and 7* are 
G-equivariant. Thus, they induce the corresponding epimorphism and monomorphism of vector 
bundles 

r ■.T*Q/G^D*, and T -.3* ^T*Q/G 

where D* := D*/G. 

The vector bundle r^, : D* ^^ Q possesses a linear almost Poisson structure. For functions 

(p,ip e C°°{D*) it is given by 

{^,V^}5. = {^or,^or}r*Q/Gor. (4.5) 

{■, ■}f), is the reduced linear almost Poisson structure that describes the reduced nonholonomic 
dynamics. 

We remark that the space of sections of the dual vector bundle r^ : D — )■ Q is naturally identified 



with the set of G-invariant vector fields on Q that are sections of D. Moreover, from (4.4) and (4.5), 
we deduce that the reduced bracket {■, -j^. is characterized by the relations 

{X,,Ye}Q, = -i7[X,Y]),, {foTQ,,Xe}f^.=Xif)oTQ„ {/ o r^., fc o r^.}^. = 0, (4.6) 

where X,Y E r(r^), and f,kE G°°{Q) while, as before, X denotes the projection of X onto Q via 



p. The main difference between the above formulae and (4.4) is the appearance of the projector CP 
in the first term. 

Local expression of the bracket {■,-}d* 

We consider a collection of vector fields {Ya} on Q that form a local basis of sections of !K on an 
open subset U O Q. Then, their horizontal lifts, {Y^}, are G-invariant vector fields on Q, that lie 
on "K, and form a basis of sections of T{t:k) on the open subset p~^{U) C Q. 

In view of the decomposition (4.3), we can complete a local basis {sj} := {Za,Y^} of sections 
of D, where the vector fields Za form a local basis of r(rvOp/G)- Note that the upper-case indices 
I, J, . . . run over the joint range of a, 6, . . . and a, (3, ... . 

A systematic way of determining Za is the following: Consider a local section s : U G Q ^ Q oi 
the principal bundle projection p : Q ^- Q. Then we may choose a basis {Za} of r(r-yi3p) along the 
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submanifold s{U), that is, Za : U C Q —t- TQ is smooth and {Za} is a basis of ^|rg)P for all q E U. 
Now, we extend the vector fields {Za} using the group action. That is, 

where (f) : G x Q ^ Q is the principal action of G on Q. Then {Za} is a basis of G-invariant vector 
fields of T^Tyop) on p~^{U). 

Next, we complete {s/} = {Za, Y^} to a local basis of sections {sj, W^} of TQ, where the vector 
fields W^ are G-invariant and form a local basis of r(r£,±). Here, D^ is the orthogonal complement 
of D with respect to the Riemannian metric S. 

Suppose that the standard Lie bracket of the vector fields s/ and sj is 

[sj, sj] = GfjSK + C}jW^. 

Then, 

T\si,Sj\=CfjSK. 

Note that the vector fields si are G-invariant which implies that Cfj and C^j are basic functions 
with respect to the principal bundle projection p : Q ^ Q = QjG. 

Let (g^) be local coordinates on Q such that 
As usual (pi) are the linear fiber coordinates on the bundle tq, : D* ^ Q, defined by pi = (s/)^. 



Then (q^pi) form a set of bundle coordinates for D*. In view of (4.6) we obtain the following local 
expressions for the bracket on D*: 

{pi,pj}d' = -Cu{q)PK, {g',Pa}5* = Kiq) {q\Pa}D* = 0' {q\f}D' = o^ 

where we have used the fact that the vector fields Y^ and Za project via p : Q ^ Q onto Y^ and 
respectively. 

Construction of the reduced kinetic energy Hamiltonian 

The restriction of the kinetic energy 9 to D induces a bundle metric 9d on D (and thus on D*). 
It is clear that this metric on D* is G-invariant and, therefore, it induces a bundle metric on D*. 
The kinetic energy associated with it defines the constrained reduced Hamiltonian function whose 
Hamiltonian vector field with respect to the bracket {■, -j^, describes the reduced dynamics of the 
nonholonomic system. 

Local expression for the constrained reduced Hamiltonian 

The constrained reduced Hamiltonian function Hfj, : D* — )■ M is given by 

H3-M,Pl) = l9''\f)PaPb+lT^mPaPp. 

where 9ij{q'') are the coefficients of the bundle metric Sd on D with respect to the G-invariant local 
basis {sj} of T^td) and {S^'^iq'')) is the inverse matrix of (9/j(^0)- Note that 

5{Za,X^) = 0, for all a and a. 

Thus, 9""('?') = 0, for all a and a. 



UNIMODULARITY AND PRESERVATION OF VOLUMES IN NONHOLONOMIC MECHANICS 



15 



The reduced equations of motion 

Using the local expressions for the bracket {•, ■}j^, it is seen that the equations of motion are given 
by (1.2). For the sake of completeness, we we give a more explicit form of the equations of motion 



using the expression of the constrained reduced Hamiltonian function 



dPa 
dt 



rPaPb+^^^PaPfS 



dq^ 



and 



dpa 
dt 



dq'- 



-s-'ci 



(4.7) 



Yl-S^'Ci^PbPi-S'^Ci 



pPiPh 



abPcPl 



nal3 /~il 



.PPPI- 



(4. 



Note that the advantage of working with the generalized nonholonomic connection is that the 
momentum variables (p/) split into two groups [pa) and (pa)? and the equations for pa are simpler. 
The above equations coincide with the universal form (3.3) but we have the additional simplification 
arising from the condition p^ = {(?SPa}j5. = 0. 

The Lagrangian version of (4.7) and ( |4.8 ) was considered in [26j (see also [IS])- If the dimension 
assumption holds, equations (4.7) and (4.8) were obtained in [23] (Theorem 3.2) as the reduced 



Hamilton equations and the momentum equation on the Hamiltonian side (for the particular case 
when the Hamiltonian function is of kinetic energy type). This Hamiltonian momentum equation 
corresponds to the nonholonomic momentum equation developed in [1] on the Lagrangian side. 



For the above reasons, equation (4.8) can be called the generalized nonholonomic momentum 
equation on the Hamiltonian side. 

The reduced metric on Q = Q/G 

We end this section by noticing that using the generalized nonholonomic connection, we may define 
a Riemannian metric Sred on the reduced space Q = Q/G. Indeed, such a metric is given by 



S{{% 



\H 



\H-s 



qjq 5 y'^'qjq )i 



for q E Q, and p{q) = q, 



(4.9) 



q, Vg G T^Q and ^ : T^Q —^ TgQ is the horizontal lift induced by the generalized nonholonomic 



^redyUq, Vg 

where U: "" ^ ^ '^ "'^'^ ^ 

connection, ^red is called the reduced Riemannian metric. Note that S^ed is well-defined since the 

metric S is G-invariant. 

This metric can be used to build a basic volume form on D 
appendix) . 



(see Remark 



4.2 



below and the 



4.2. Unimodularity of the reduced linear almost Poisson structure. The constructions in 
the previous section will be useful to study the unimodularity of the reduced almost Poisson structure 
{•, ■}f,,. Recall from (4.2) that the reduced space D = D/G can be decomposed as 

3 = Vp/G © :k/g. 

Our study of the unimodularity of the almost Poisson structure {■, -jg. will be performed by 
separately considering linear functions induced by sections of the spaces involved in the above de- 
composition. 
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Assume that both the manifold Q and the vector bundle D* are orientable. Then we can choose 
a basic volume form $ on D*. 

Remark 4.2. The volume form $ can be chosen in different ways. For instance, we can take the 
Riemannian volume form i/g^^^ on Q associated with the reduced metric Sred and the volume form 
A^^ on the fibers of D associated with the bundle metric 9d on D which is induced by the restriction 
Sd of the G-invariant metric S to D. Then, following the notation in the Appendix, $ = z/g^^^ A A^^ 
is a basic volume form on D*. o 

Denote by M$ the modular vector field of D* with respect to the basic volume form $. We define 
the family "Ji of 1-forms Qq, on Q that satisfy 

w$(r)or5, =M$((F%), (4.10) 

where Y is any section of "K and {Y^)(, den otes the linear function on D* induced by the section 



Y^ G r(r^). Note that by Proposition 3.3, the vector field M<j, is rg*-vertical. This implies that 
M$((y^)£) is a basic function with respect to Tj^^ and the relation (4.10) is consistent. 



If the dimension assumption holds, then CK = TQ and a}$ is uniquely defined. In other cases, 'Jl is 
the space of sections of an afiine subbundle of TQ with rank 

dimQ — rank(D + Vp). 

We are now ready to give necessary and sufficient conditions for the unimodularity of D in terms of 
the family "R of 1-forms w$. 

Theorem 4.3. Let $ he any basic volume form on D* . The linear almost Poisson structure {■, -jj^* 
on D* is unimodular if and only if the following conditions hold: 

(i) For every section Z ofV^p/G we have 

div$(XzJ = 0, 

where Xz^ is the Hamiltonian vector field on D* associated with the linear function Zi : D* ^^ 
M induced by the section Z . 



(ii) There exists a 1-form a}$ on Q belonging to the family Jl defined by (4.10) such that 

a}$ = —da, 
for a certain smooth function a on Q. 



3.6 



Pro of. Suppose that the almost Poisson structure is unimodular. From Remark |3.5| and Theorem 
it follows that there exists a G C^IQ) such that 



M^ = X^,^^^. (4.11) 

Let Z as in item (i) of the statement of the theorem. Then Z is a G-invariant vector field on Q that 



lies on the intersection V^p = VpCi D. Using the definition of the bracket on D* given in (4.6) we 
find 

X^oTgAZe) = {Ze,a o t^,}^, = -Z{a) o r^.. 
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where Z is the projection (via p) of the vector field Z. Thus, since Z is p-vertical, we have 
XsoT-X'^t) — 0- Therefore, as a consequence of (4.11) we obtain that the condition (i) in the 
theorem is necessary for unimodularity. 

Similarly, if F is a vector field on Q that belongs to the distribution 'K then, since the 1-forms in 



'51 satisfy (4.10), from (4.11) we obtain 

u:^{Y) o r^. = M$((F^),) = X^o,^, ((F^),) = {(F^),, a o r^ J 



D*-l D* 



-r(a)or5. 



where we have used again the definition of the bracket on D* given in (4.6). The last equality 
implies that a}$(y) = —Y{a) = —da{Y). This shows that item (ii) in the theorem is also a necessary 
condition for unimodularity. The sufficiency of these conditions is a consequence of the decomposition 
dO). D 



Local expressions for unimodularity 

We are interested in obtaining a local version of the above theorem to apply it to concrete examples. 
We recall the notation that was used above: {q'') are local coordinates on Q, whereas {pa,Pa) are 
linear coordinates on the fibers of D* induced by the sections Za, Y^ that are a basis of D adapted 
to the decomposition (4.3): 

D = v^p/G © :k/g. 

Consider the basic volume form ^ = u A A on D* and assume that with our choice of coordinates 
we have 



;\dq\ A = /\Z„/\f/. 



Then 



$ = yy dq' /\ dpa /\ dpa . 



(4.12) 



On the other hand, using ( |3.2[ ) and the fact that {Za)e = Pa and {Y^)i = pa we get the following 
expressions for their Hamiltonian vector fields 



*-■ = icLp.)^^HcLP,)^^ 






^-- = yi^t + KaPi)i- + K,p^) ^ 



'dq' 



dpa 



dpp 



(4.13) 
(4.14) 



Thus, from (4.12) and (4.13), we immediately get the following local version of item (z) in Theorem 

(4.15) 



Ci, = Cl + C. = for all a. 



Similarly, from (410), (412) and (414), it follows that 



Uq,{Ya) 



dYi 



+ ^aa + ^afz 



(416) 
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Therefore, we conclude that condition (ii) in Theorem 4.3 is equivalent to the following equations 
for the unknown a G C°°{Q): 



-Y„ 



o] 






+ ^'. + a 



'/3 
a/3 



for all a. 



(4.17) 



4.3. Invariant measures for the reduced nonholonomic system. We now state our results 
looking ahead at measure preservation for nonholonomic systems. 

Corollary 4.4. The red uced nonholonomic dynamics preserve a volume form if and only if conditions 
(i) and (ii) in Theorem 4-3 hold for a basic volume form $ on D* and a real C°° -function a on Q. 



Furthermore, if these conditions hold then the dynamics preserves the basic volume form $' on D* 
given by $' = e^°'^s»(|). 



Proof. The first par t follows from Theorems 3.8 and 4.3 On the other hand, if conditions (i) and 
(ii) in Theorem 4.3 hold then the modular vector field M$ in D* is just the Hamiltonian vector field 
of the basic function a o tq^,. Thus, using (2.4), we deduce that M$/ = 0, which proves the second 
part. n 

We will now specialize our discussion to different reduction scenarios. This will allow us to express 
conditions (z) and {ii) in Theorem 4.3 in simpler terms that are useful in applications. In doing this 



we shall recover several results existing in the literature in a unified manner [HI [20], |25l H2] . 
Invariant measures under the dimension assumption 



Suppose that dimension assumption (4.1) holds, that is, 

Vp + D = TQ. 

Then the horizontal space S) of the nonholonomic connection equals "K, and the distribution [K = TQ. 
As a direct consequence of Corollary 4.4 we recover the characterization of systems with an invariant 
measure given in |12] in a local form. Here, we present the intrinsic version: 



Corollary 4.5. [12] The reduced nonholonomic dynamics preserves a volume form if and only if 
there exists a basic volume form $ on D* such that the following conditions hold: 

(i) For every section Z ofV^p/G we have 

div$(Xz,) = 0. 



(ii) The 1-form cD$ on Q defined by (4.10) is exact. 
Furthermore, if (i) and (ii) hold, and uj<^ = —da, then the dynamics preserves the volume form 



$' 



''d*<|). 



The local version of condition (i) in the above corollary is again given by (4.15). However, the 



local version of condition (ii) given in (4.17) greatly simplifies in this case. Since CK = TQ, the local 

d 



vector fields Y^ can be chosen as 



Y. 



(9g"' 
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where (g") are local coordinates for Q. Equation (4.16) implies that w$ is locally given by 



C^=(ci + Cf^)rfr- (4.18) 



Invariant measures for G-Chaplygin systems 

Now, suppose that the nonholonomic system is a Chaplygin system, see e.g. [2JJ, that is, 

TQ = Vp(B D. 

In this case the distribution V^p = and the condition (i) in Theorem 4.3 is empty. On the other 
hand, condition (ii) yields the following result that was proved in |8]. 

Corollary 4.6. ^ For a G- Chaplygin system, the reduced nonholonomic dynamics preserves a vol- 
ume form if and only if the 1-form a)$ on Q is exact. Furthermore, if u)<s> = —da, with a e C°°{Q), 
then the dynamics preserves the volume form 



In this case, the local expression (4.18) for a)$ further simplifies to give 

oj<, = C^^l,dr. (4.19) 

Invariant measures for LL systems 

Now suppose that the nonholonomic system is LL. This means that the configuration manifold 
coincides with the symmetry Lie group, i.e. Q = G. In this case, both the metric S and the constraint 
distribution D are left-invariant. Hence they are respectively determined at the Lie algebra g by a 
symmetric, positive definite, inertia tensor J : g — )■ g*, and a subspace 9 C 3 that is not a Lie 
subalgebra. For all /i G G we have 

9{h){uh,VH) = (l((T,/,-i)M),TA-i(t;h)), D{h) = {T,k){d), 

where Uh, Vh G T^G, e is the identity element in G, (-, ■) denotes the duality pairing, and the mapping 
//i : G — 7- G is left multiplication by h. 



For an LL system, the condition (ii) in Theorem 4.3 is empty. On the other hand, the local 



expression 4.15 of condition (i) simplifies to: 



Cl^ = 0, for all a. 

Moreover, in this case the vector fields Za can be taken as left extensions of a basis {ca} of d (i.e. 
Za{h) = Tclh{ea) for h E G). Then if we denote by IP : — )■ 5 the orthogonal projector corresponding 
to the metric defined in q by the inertia matrix J, the following result follows 

Corollary 4.7. The reduced nonholonomic dynamics preserves a volume form if and only if 

n 

^(e^y[e„e4) = 0, for all a, (4.20) 



where {ca} is a basis of d and {e"} is the dual basis of d* . In addition, if (4. 20) holds then the 
nonholonomic dynamics preserves the euclidean volume on d* . 

The conclusion about the euclidean volume is established by noticing that any volume of basic 
type on d* is a constant multiple of it. 
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Remark 4.8. If the co-dimension of (9 in g equals to 1, then the condition (4.20) can be given in 



simpler terms. Let /3 G 0* span the annihilator 9° of d, and let rj := I ^/3 G 0. Then, it is easy to 



prove that condition (4.20) is equivalent to the requirement that 



— -— ad*/3 + Mg = yU/3, 

for a certain yU G M, where ad* : g x g* — )■ g* is the coadjoint representation of q, and Mg G Q* is the 
modular character of g (we recall that Mg(,^) = trace(ad5), for ^ G g) . This result had been proven 
by Jovanovic [20] (see also |25]). o 

5. Algorithm to investigate the existence of invariant measures 
In this section we give an algorithm to determine if a nonholonomic system with symmetry pos- 



sesses an invariant measure. The algorithm is written so that the conditions in Theorem 43 can be 
applied in practice in a systematic way. We restrict our attention to a more specific kind of systems 
that satisfy the following conditions: 

CI. The first de Rham cohomology group of the shape space Q = Q/G is trivial. 

C2. There exists an open dense set U G Q with a global chart that defines coordinates (g^). 

For the examples that we treat in the next section, the space Q is M^, S"^ or 5*0(3), that satisfy 
both conditions CI and C2. We recall that in our developments we are always assuming that the 
intersection V^p := D nVp has constant rank. 

The steps are described above under the assumption that (S,-D) is a symmetric nonholonomic 
system on Q. Most of the steps only involve computations that can be systematically performed 
with a symbolic algebra program. 

Step 1. Find a basis {Za, Y^} of G-invariant vector fields of D in such a way that {Za} is a basis of 
r(rvcp) and {Ya} is a basis of T{t^). 

The vector fields Za and Y^ need not be defined globally. It is sufficient that they are defined 
on the dense open set p~^{U) C Q. In the examples treated below, the vector fields {Za, Y^} 
were found by inspection. In Section 4.1[ we presented a systematic way of determining the 



vector fields Za and Y^. 

If the dimension of V^p is zero, skip steps 2 and 3 and go to step 4. On the other hand, if 
the dimension of CK is zero, the algorithm terminates in step 3. 
Step 2. Compute the structure coefficients Cai defined by 

7 {[Za, Z,]) = CiZ, + C:,Y^, 7 {[Za, rj) = ClZ, + C^aaY^, 

where CP is the S-orthogonal projection onto D and [-, ■] is the standard commutator of vector 
fields. Notice that by G-invariance of the basis {Za,Ya} and the metric S, the structure 
coefficients Caj sue functions of g'. 
Step 3. A necessary condition for the existence of an invariant measure (coming from item (i) in 
Theorem 4.3 and expressed locally in (4.15[)) is that 



Cat + C. = 0, for all a. 



If the dimension of "K is zero, this condition is also sufficient. 
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Step 4. Compute the coefficients C^^ defined by 

y{[Y.,Yp\)=Cl,Z, + ClpY,. 

As in step 2, here one needs to compute the standard commutator of vector fields and their 
9-orthogonal projection onto D. As before, one should interpret the coefficients C^j as 
functions of g^ 
Step 5. Consider the projections Y^ (via p) on U of the vector fields Y^, and write them as 

Y = P— 
Now, denote by "R the family of 1-forms a) on [/ that satisfy 



dY^ 



A necessary condition for the existence of an invariant measure is that there exists a 1-form 
a) in the family 3? that satisfies uj = —da^ for a certain dfj G C°°{U). (In particular notice 
that such 1-form cD satisfies du = 0.) 
Step 6 Let {q^PaiVa) be a system of local coordinates on D* induced by the basis {Za,Ya}. Then 
the volume form 

$ = e^e l\dq' l\dpa l\dpa 

L a a 

is preserved by the system on the open set {'rg^)^^{U) C D*. One should finally verify that 
the above volume form admits an invariant extension to D*. 

Remark 5.1. If the dimension assumption holds, then such an invariant extension always 
exists. To see this fix a basic volume form \1/ on D*, and suppose that 

vl>|(^^^).,(^) = e^^oorn^^ /\dq' /\dpa /\dp^, with (^ G C^{U). (5.1) 

t a a 

In such a case, there is a smooth extension ^ G C°°{Q) of the real function a^ — (fj and 

v[/' = e(^°^5*)vl> 
is an invariant volume form. 



In fact, if cDij, is the (global) 1-form on Q given by (4.10) then a direct computation using 



(5.1) proves that 

{Q^)\Q = + dCQ = -d(aQ-(^). (5.2) 

In particular, w^ is closed in U and, since U is dense, we obtain that w* is closed (in Q). 
This, using CI, implies that oj<!j is exact, that shows that condition (ii) of Theorem 4.3 is 



satisfied. o 
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6. Examples 

6.1. Body with planar section rolling over a fixed sphere. Following |27l EB], consider the 
motion of a rigid body that possesses a planar facqj that rolls without slipping over a fixed sphere 
of radius R. We consider a body frame {£"1, £"2, £3}, whose origin is located at the center of mass C 
of the body, and such that the E's-axis is normal to the planar face of the body. The spatial frame 
{ei, 62, 63} has its origin at the center O of the sphere. The distance between the center of mass C 
and the planar face will be denoted by L We assume that our choice of body frame is such that E^ 



is the outward normal vector to the sphere at the contact point P (see Figure 6.1). 

The two frames are related by an attitude matrix g G 5*0(3). Let a; G M'^ be the spatial coordinates 

of the vector OC . The body coordinates of OC are X = g~^x. At any configuration, the condition 
that the planar face of the body is tangent to the sphere is expressed by the condition 

X3 = R + i. (6.1) 




Figure 6.1. Body with planar section rolling over sphere 

The constraint that the planar face rolls without slipping is equivalent to the requirement that the 
contact point P is at rest. This is expressed by the relations 

Xi = -Rn2, X2 = i?fii, (6.2) 

where the time derivative of the coordinates Xi,X2 is computed in the body frame and the vector 
n = {Qi, Q2, ^3) £ ^^ is the angular velocity of the body expressed in the body frame. We think of 
f2 as an element in the Lie algebra so(3) corresponding to the skew-symmetric matrix g^^g with the 
usual identification of so(3) with M^ via the hat map (see e.g. pil]). 

In view of our discussion, we can take Q = S0{3) x M^ with coordinates (Xi,X2) for M^ as the 
configuration space. Indeed, if we know the attitude matrix g, and the value of the coordinates 



an everyday life example of such type of rigid body is a shoe (without heel). 
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(Xi,X2), then, in view of (6.1), we can specify the configuration of the body. The point of contact 
P is easily determined since the vector OP has space coordinates Rg ( ° ) • 

Let Ci , 62 , 63 be the left invariant moving frame of S0{3) obtained by left translation of the 
canonical basis at the group identity. We have the commutation relationships: 

The entries fli, Q2, ^3 of the angular velocity in the body frame fl are quasi-velocities with respect 
to this moving frame. 



The no-slip constraints (6.2) define the constraint distribution D C TQ that is generated by the 
vector fields 

The kinetic energy of the body is given by 

1 777 

x = -{m,n) + -\\x\\\ 

where the derivative x is computed on the space frame, (■, ■) is the euclidean scalar product in M^, 
m is the mass of the body, and the 3x3 symmetric, positive definite matrix I with entries Jjj is the 
inertia tensor of the body. Notice that in view of our choice of body axes, I need not be diagonal. 
However, by an appropriate rotation of the axes Ei, E2, we can assume that /12 = 0. 

Using that g G 5*0(3) we have 



I . I i2 II — 1 • I i2 



and hence, since 



we can write 



g-^x = X + flx X, 

1 777 / • \ 777 

x = -{in,n) + —\\n X xW^ + m(x,n X X) + —wxW^. (6.4) 



Denote by {e^,e^,e^} the dual basis to {61,62,63}. A straightforward calculation using (6.1) 



shows that the kinetic energy % defines the following metric S in Q, 



)e^ ® e^ 



g = {hi + m{Xl + {R + if))e^ ®e^ + {h2 + m{Xl + {R + if))i 



+ (/33 + m{Xl + Xl))e^ ® e^ - 2mXiX2e^ ® e^ + 2(/i3 - mXi{R + £))e^ ® e^ 

^ ^ ^ ^ (6-5) 

+ 2(/23 - mX2{R + i))e^ ®e^+ 2m{R + i)e^ ^ dXi - 2mX2e^ ® dXi 

+ 2mXie^ (g) dX2 - 2m{R + i)e^ ^ dX2 + m{dXi (g) dXi + dX2 6^X2). 

The group G = 50(3) acts on Q by left multiplication on the 50(3) factor and leaves the metric 
S and the constraint distribution invariant. The space V^p has dimension 1 and is generated by 63. 
The shape space Q = Q/G = M^ that has a global chart and its first de Rham cohomology group is 
zero. Hence, we are in the framework to apply the algorithm described in section |5| 
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The vector fields 

Za=l = 63, Ya=i = - — 62 + -=—1 ■ ,^^0 , ^^n,, 63 + 



R \R{l33 + m{Xf + Xl))J dXi 

-In + mix, \. d 



"=' R' Whs + m{Xf + XI)) J ^' + 



9X2' 
satisfy tlie properties of step 1 of the algorithm. To perform step 2, we compute the commutators: 

[Za=l , Ya=l] = -^Ci , [Za=l , Ya=2] = "^62. 

Next, we need to calculate the S-orthogonal projection of these vectors onto D and express them 
in terms of the basis {Za=i , ^0=1 , ^0=2}- This is done by solving the linear system of equations 
TVa = ba, where 

/ g(z,=i, Za=i) \ / g(z„=i, [Za=i , rj) 

T=\ g(r,=i,r«=i) g(r„=i,r„=2) , h^=\ g(ra=i,[Za=i, >^J) 

The entries y^ of ija are the structure constants Cf^]^„. After solving the systems of equations 
corresponding to a = 1, 2, with the aid of MAPLE"'"'^, we obtain: 

777- 

C:=l^^, + C=i'«=2 = ^3det(T) (-^23(/ii + me)X, + /i3(/22 + me)X2 + m£(/n - /22)XiX2) . 

According to step 3 of the algorithm, a necessary condition for the existence of an invariant measure 
is that the above expression equals zero. Since such equality must hold for all values of Xi and X2, 
we obtain the necessary conditions for the existence of an invariant measure: 

(Jn + m£2)/23 = 0, (J22 + mi-')h3 = 0, {hi - 122)^ = 0. 

The first two conditions imply that /13 = /23 = (i.e., the inertia tensor is diagonal), since Jn, /22 > 
0. The third condition is satisfied if either £ = or In = /22- 

The above conditions are also sufficient for the existence of an invariant volume. Indeed, the 
measure for the first case was found in |41J and for the second case already in [37] (see also [B]), 
where the equations of motion were also explicitly integrated. 

Therefore we have: 

Theorem 6.1. The reduced equations for a rigid body with a planar face that rolls without slipping 
on a fixed sphere possess an invariant measure if and only if at least one of the following conditions 
hold: 

(i) the center of masses is in the base plane (£ = 0), which actually implies that the body is flat, 
(ii) the body is axially symmetric (In = I22, hs = I23 = 0). 

The preserved measures have the form 

(/n/22 + mluX^ + ml22Xl) ■ {h-i + m{Xl + Xl))^'^ dXi A dX2 A dQi A rf^s A dQs, 
and, respectively, 

(In + m{X^ + Xl) + mf) ■ {hiihs + m{Xl + X^)) + mfh^)^/'^ dX^ A dX2 A d^i A d^s A ^^3. 
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6.2. The Chaplygin top. Consider the motion of an inhomogeneous ball that rolls without slipping 
on the plane. If the center of mass of the sphere coincides with the geometric center, we recover 
the well known problem of the Chaplygin sphere, that possesses an invariant measure [10]. Here, we 
consider the general case that is sometimes referred to as the Chaplygin top |31j. In the case of an 
axisymmetric ball, this problem was originally studied in detail in |j32j and [TT]. 

Let R be the radius of the sphere. The space frame {ei, 62, 63} is chosen so that the rolling takes 
place on the plane z = —R. We consider a body frame {i?i, £"2, £^3}, whose origin is located at 
the center of mass of the sphere M, and such that the geometric center C of the sphere lies on the 
£'3-axis. Let < £ < Rhe the distance between the geometric center of the sphere and the center of 
mass. The body coordinates of the center C of the sphere are (0, 0, —i). See Figure 6.2 below. 




Figure 6.2. Chaplygin Top 

The configuration space of the system is Q = SO (3) x M^. A matrix g G SO (3) specifies the 
orientation of the ball by relating the body and the space frame. We will use Euler angles as local 
coordinates for 5*0(3). We use the x-convention, see e.g. [30] and write a matrix g G 5*0(3) as 

(cos ip cos (f — cos 9 sin ip sin ip — sin ^ cos ip — cos 9 sin (f cos ^ sin 9 sin (f 
cos ip sin if + cos 9 cos ip sin ip — sin ^ sin (p + cos 9 cos (p cos ^ — sin 9 cos ip 
sin 9 sin ip sin 9 cos ip cos 9 

where the Euler angles < ifyip < 27r, < ^ < vr. According to this convention, we obtain the 
following expressions for the angular velocity in space coordinates uj, and in body coordinates f2: 

(9 cos ip + tp sin ip sin ^ \ / 9 cos ip + tp sin ip sin 9 \ 

9 simp — Ip cos ipsm9 , f2 = I — ^ sin ■?/' + cos ^/^ sin 6' . (6.6) 

ip + Ip cos 9 / \ ip cos 9 + ijj ) 

Let (x, y, 0) be the spatial coordinates of the geometric center C of the sphere. The constraints of 
rolling without slipping are 

i, = Rl>J2 = R{9 svmp — il) cos ipsvuB)^ y = —RxjJi = —R{9 cos ip + ip simp sin 9). (6.7) 

The kinetic energy of the sphere is given by 



1 m 



u 



(6.8) 



where u are the space coordinates of the center of mass. Here, I is the inertia tensor of the sphere 
with respect to the center of mass with entries lij i,j = 1, 2, 3. Since we have already made a choice 
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of the body frame, we cannot assume that I is diagonal. However, by an appropriate rotation of the 
body frame around the E's-axis, we can assume that /12 = 0. 

The expression for u in our coordinates is given by 

(X \ / \ / X + £ sin 6* sin (y9 

y \+g\ = y-isiwOcos^ \ . (6.9) 

y \ £ / \ loose 

Therefore 

(X + iO cos 9 sin (p + lip sin 9 cos (p 
y — £9 cos 9 cos ip + iip sin 9 sin (p 
-i9 sin 9 

The local expression for the kinetic energy metric S of the problem is obtained by substituting the 
expressions for u and f2 onto (6.8). 

Symmetries 

There is a freedom in the choice of origin and orientation of the space axes {61,62}. This corre- 
sponds to a symmetry of the system defined by a left action of the Euclidean group SE{2) on Q. 
Let 

(cos tD — sin ^ V 
sin d cos d w 
1 

denote a generic element on SE(2). The action of /i on a point q E Q with local coordinates 
{ip,9,'i/j,x,y) is given by 

h ■ q = {{p + {},9,iIj, xcos-i? — y sm{} + v, x sm{} + ycos-d + w). 

One can check that both the constraints and the kinetic energy are invariant under the lift of the 
action to TQ. The action of SE{2) on Q is free and proper and the shape space Q = Q/G = S"^. In 
our local coordinates the orbit projection p : Q ^ S'^ is given by 

p{ip,ilj,9,x,y) = {9,^/j), 

where {9, ijj) are spherical coordinates on the unit sphere 7i + 7! + 73 = 1; defined by 

7i = sin6'sin'?/', 72 = sin^cos-?/', 73 = cos 6'. 

The vector 7 = (71,72,73) = 5'~^63 is the the usual Poisson vector whose entries are the body 
coordinates of the unit vector that is normal to the plane of rolling, and evolves according to the 
kinematic equation 

7 = 7 X f2. 
The vertical subbundle Vp is spanned by 

Odd 



Vp = span , , 

op> ox oy 
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On the other hand, the constraint distribution D is spanned by the SE(2)- invariant vector fields 

Z - ^ 

d ^ d ^ d 

Xa=i = 7^ + Rsmip- Rcosif—- , 

oti ox ay 

d d d 

Xa=2 = TTT — -R COS (y9 sin 6^7- RsmifismO-—. 

dip ox Oy 

It is then clear that the intersection V^p = D (1 Vp has constant dimension 1 and is spanned by 
Za=i- Moreover, since S"^ is simply connected, its first de Rham cohomology group is zero. Also, the 
angles (6','?/') form a coordinate chart on an open dense set of S"^, so all of the conditions to apply 
the algorithm in Section |5] hold. 

The following vector fields, together with Za=i satisfy the requirements of step 1 of the algorithm: 

Ya ■— y^a — 777^ ^ T ^a=\ a — L, Z. 

b[Za=l, Za=l) 

We compute the standard commutators: 

d d 

[Za=i,Ya=i] = Rcosip-- + RsmLp—- + XiZa=l , 

ox oy 

d d 

[Za=i,Ya=2] = -Rsinv9sin6'- i?cosv9sin6'-- + A2^a=i, 

ox oy 

where Ai,A2 are certain functions of {0,ip). We should now compute the projection of the above 
vector fields onto D and express them as a linear combination of Za=i,Ya=i,Ya=2 to determine the 
coefficients Cf^^ ^. In fact, looking ahead at step 3 of the algorithm, we are interested in computing 
C^^i „ for a = 1, 2. A simple linear algebra argument shows that C^^i^^^ coincides with the 
component of Xa=i when the orthogonal projection of 

9 J, . d 
Rcosip- — h Rsimp—- 
ox oy 

onto D is expressed in terms of the basis Za=i, Xa=i, Xa=2- A similar idea can be used to compute 
Ca=ia=2- Using these observations and with the aid of MAPLE"""^ we obtain: 

where T is the (positive definite) matrix 

S(-^a=l,^a=l) S(^a=l,-^a=l) S(^a=l, -^0=2) 
S(^a=l,-^a=l) S(-^Q=l,-^Q=l) S (-^0=1,-^0=2) 
5{Za=l,Xa=2) S (-^0=1,-^0=2) 9{Xa=2,Xa=2] 
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and the coefficient functions ai, 02, &i, &2 are given by: 

am = -/i3 r^^ + (/22 + m{R' + f)) cos(^) + ^ cos(20) 

02(6') = /i3/23sin6', 

hm = /23 (^ + (/n + miR' + f)) cos(^) + ^ cos(20)) , 

HO) = ^ ((^22 - /ii)(/33 + mR'' + miRcosO) - 4 + 4) . 
The necessary condition for the existence of an invariant measure, coming from step 3 of the 



algorithm, is that the expression (6.10) vanishes identically for all {9,ip) in the chart, that is, for 
all < ^ < TT, < '?/' < 27r. By linear independence of cos-?/', cos(2'?/'),sin'?/^,sin(2'?/'), it follows that 
this can only happen if either i = or all of the coefficient functions 01,02, 61,^2 vanish. A quick 
examination of the above expressions, shows that the latter case can only hold if /23 = /13 = and 

-^11 = -^22- 

That these conditions are also sufficient is proven in [JT] (see also [6j). There it is shown that the 
reduced equations of motion can be presented in the vectorial form 

K = Kxn + mR{j xfl) X {fix p), j = j x ft, (6.11) 

and possess the invariant measure: 

,, , , ^= = {hi + m||p|n V//11/33 + m{p, Ip) dj A dn 

y/Iuh3 + m{p, Ip) 

In the above formulae p is the vector that connects the contact point with the center of mass of the 
sphere written in body coordinates, and K is the angular momentum of the ball with respect to the 
contact point, also written with respect to the body frame. Explicitly we have: 

p = Rj+\ , K = in + mpx {fix p). 

Therefore we have: 

Theorem 6.2. The reduced equations for the Chaplygin top possess an invariant measure if and only 
if at least one of the following two conditions hold: 

(i) The center of mass of the sphere coincides with the geometric center (i = 0). 
(ii) The ball is axially symmetric (In = 122,-^13 = -^23 = 0). 

We stress that the above conditions were known to be sufficient but we have just shown that they 
are also necessary. 

6.3. The rubber axially symmetric Chaplygin top. We will now show that an invariant measure 
also exists when one adds a rubber constraint to the symmetric Chaplygin top and we will use our 
method to provide an explicit formula for such an invariant measure. The rubber constraint prohibits 
the ball to rotate about the axis 63 that is perpendicular to the plane where the rolling takes place 



see Figure 6.2). This same problem, for a non-symmetric ball, is considered in [5] where the authors 



show non-existence of an invariant measure by exhibiting the presence of limit cycles in the dynamics. 
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We will keep the notation of the previous section and work under the assumption that the ball is 



axially symmetric (In = 122,-^13 = hs = 0). The rolling constraints (6.7) are complemented by the 



rubber constraint (see (6.6)): 

Ws = + tp cos9 = 0. 

Symmetries 

The same G = SE{2) symmetry exists with vertical space 

d d d 



Vp = span ^ . . 

Oif ox oy J 

However, in this case the constraint distribution D is spanned by the G-invariant vector fields 

d d d d 

Ya=i = — cos6- — h 7— — i?cosv9sin^7- RsinipsmO^— , 

dip dip ox oy 

9 j^ . d ^ d 

>a=2 = 7^ + -R Sm V?7^ -K COS if—- . 

OO ox oy 

Hence we have D fl Vp = {0}. Moreover, we have TgQ = Dq(B Vpq for all g G Q so we are dealing 
with a Chaplygin system. 

The vector fields Ya=i,Ya=2 defined above verify the requirements of step 1 in the algorithm. 
Moreover, the condition in step 3 is vacuously satisfied. We thus concentrate in examining the 
condition in step 5. 

The projections Y^ = Tp{Ya) are given by 

Y -± i> -A 

We consider the 1-form Q that is defined by the relations 

We compute the standard commutator 

d 

[Ya=i,Ya=2] = -sin6' — . 

Using MAPLE^'^ we find that its orthogonal projection onto the constraint distribution D is 
f{9)Ya=i where 

, , _ sin g((/ii + me'^ - 133) cos 6 + mJR) 

■'^ ' ~ (111 + m£2) cos2 e + hs sin^ 6 + 2miR cos 6 + mR^ ' 
It follows that 

Therefore Q = —f{6)dd = —da where 

a{9) = --In {{hi + mf) cos'^ 9 + hssin^ 9 + 2miR cose + mR^) . (6.12) 

Since S"^ is simply connected and the spherical coordinates ip, 9 form a chart on an open dense set, 
it follows that an invariant measure exists. Below we explain how to express this measure in terms 
of physical variables. 
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Expression for the invariant measure in physical variables 

We will write the invariant measure in terms of the angular velocity 17 written in body coordinates 
and the Poisson vector 7. 

We start by considering the vector fields 

^ TT. ^ • ^ TT. . d d 

PVi := 7—, Py2 := cosv^T^ + smv9--, PVs :=- sm </)-- + cos v^tT' 

0(f ox oy ox oy 

that form a G- invariant basis of a complement of D in TQ. 

The vector fields Ya=i, Ya=2, Wi, W2, W3 form a local basis of sections of TQ. Denote by f 1, . . . , fs, 
the quasi-velocities associated to this basis. We have: 

vi='ip, V2 = 0, V3 = i^ + ip cos9, V4^ = Rip sm9 + x cos ip + y simp, v^ = R9 — x simp + y cos ip. 

The constraint distribution D is described by the relations f 3 = f 4 = f 5 = while f 1, f 2 are natural 
coordinates for its fibers. 



Now we need to pass to momentum variables. When the kinetic energy (6.8) is restricted to D, it 
can be written in terms of the quasi-velocities vi and f 2 as 

X =- ((sin^ eiihi + mf) cos^ 6 + /gg sin^ 6 + 2miRcos 6 + mR^)vl 

+ ((/ii + m{R^ + £2) + 2mRicose)vl) . 
This yields the following expressions for the quasi-momenta pi,P2'- 

Pi = (sin^ 0{{lii + mf) cos^ 6 + I33 sin^ 6 + 2ml R cos 9 + mR^)vi, 
P2 = (111 + m(i?2 + f) + 2mRicos9)v2. 

According to step 6 of the algorithm, the local expression for the preserved measure is $ = 
e'^dip A dO A dpi A dp2- In view of (6.12) this is 

$ = ((/^^ + mf) cos^ 6 + I33 sin^ + 2miR cos 6 + mR^) ~^^^ dip AdO A dpi A dp2 
= sin^ 9{Iu + m{R^ + f) + 2mR£ cos 6) 

■ {{hi + mf ) cos^ 6 + /33 sin^ 6 + 2m£i? cos 9 + raR^) ^''^ dip A d9 A dvi A dv2. 

Below we shall consider the equations of motion in vectorial form. In such a formulation, the 
vertical component of the angular velocity u^, = (7, fl) = ^3 is a constant of motion. Therefore we 
append the equation ^3 = to the system, that will now preserve the measure 

!> = sin^ 9{Iii + m{R^ + (P^) + 2mR(i cos 9) 

■ {{hi + mf ) cos^ 9 + J33 sin^ 9 + 2m£i? cos 9 + raR^) ^''^ dip A d9 A dvi A dv2 A dv^. 



Now notice that, in view of (6.6), the components of ft are related to fi,f2,f3 by 

- sin ip cos 9 sin 9 cos ip sin ip sin 9 

- cos Ip COS 9 sin 9 — sin ip cos ip sin 9 
sin^ 9 cos 9 
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The determinant of the above matrix is sin 9 and thus 
$ = sin^(Jii + m{R^ + f) + 2mm cos 9) 

■ ((/ii + mf) cos^ 9 + /33 sin^ 9 + 2mm cos 9 + mR^) ^'^ dip Ad9 A dQi A dQi A dQs. 

Finally notice that sin 9d'ipAd9 is the area form on the sphere 1 17| P = 1 induced from the Euclidean 
metric d'-yi A ^72 A ^73 on M^. Using these considerations, and putting cos 6^ = 73, and sin^ 6^ = 7^ + 7I 
we define the following volume form on M^ 

y, = (111 + m{R^ + f) + 2mRi-f3) 

1/9 (6-13) 

■ ((/n + mf)^l + 133(72 + 72) + 2m^i?73 + mR^ d-y A dn. 

The vector form of the equations that we present below preserves the above measure when restricted 
to the sphere II7IP = 1. 

The equations of motion in vector form. We now write the equations of motion for the symmet- 
ric, rubber Chaplygin top in vector form. According to the Lagrange-D'Alembert principle, when 



the rubber constraint is enforced, the equations of motion (6.14) are modified to give 

K = K X ri + mi?(7 X ri) X (ri X p) + A7, 7 = 7x^2, (6.14) 

where the multiplier A is uniquely determined from the condition uo^ = (^2,7) = 0. Differentiating 



this relation one obtains (fi,7) = 0, that combined with (6.14) yields after some algebra 



{Kxn + mR{j xn) X {fix p)- 2mR{p, 7 x n)n + mR{p, n)j X ft , I{p)-^j ) 

~ a(p)-^7,7) 

where J(p) denotes the 3x3 matrix 

J(p) = I + m| IpI |2 ■ ^^3x3 - mpp^ 
that relates K and fi (K = J(p)ri). 



Notice that equations (6.14) have conserved quantities (r2,7) and ||7||2. Our interest is in the 
level set {fl,j) = and ||7il^= 1- We have verified that these equations preserve the measure (6.13) 
on the level set ||7||2 = 1 using MAPLE 



TM 



6.4. Dynamically balanced sphere rolling on the exterior of a cylinder. We consider the 
motion of an inhomogeneous sphere, whose center of mass coincides with its geometric center C, 
that rolls without slipping on an infinitely long vertical circular cylinder of radius r. If the ball is 
homogeneous the problem is integrable for general shapes of cylinders as was shown in [6] . We also 
mention that the problem of a homogeneous ball rolling on circular cylinder in the presence of gravity 
was considered by Routh. 

The space frame is chosen so that the z-axis coincides with the axis of the cylinder. Let i} be the 
polar angle of the center of the sphere C on the xy plane. Then, the space coordinates of C are 



{{R + r) COS"*?, {R + r) sin??, z) where R is the radius of the ball. See Figure 6.3 below. 

The configuration of the ball is completely determined by the position of its center C and by a 
rotation matrix g G 5*0(3) that relates the space frame with a fixed body frame, whose origin lies 
on the center of the sphere C and that we assume to be aligned with the principal axes of inertia of 
the sphere. Hence the configuration space for the problem is Q = S0{3) xM. x S^. Our coordinates 

for Mx S^ are (z,^). 
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Figure 6.3. Dynamically balanced ball rolling on a cylinder 

The constraint of rolling without slipping is obtained by requiring that the contact point of the 
ball with the cylinder is at rest. This gives 

'& = { — 1^3, i = -Rfwi sin ^ — W2 cost?), (6.15) 

\R + rJ 

where u: = {ui, U2, Ws) is the angular velocity of the sphere written in space coordinates. In terms of 
the Euler angle convention introduced in the previous examples, we have 

R 



^ 



R + r 



if + cos 9ip 



z = R( — sin((y9 — '{})9 + sin 9 cos{ip — '{})ip ) . 



Since the center of mass of the sphere coincides with the geometric center, the kinetic energy is 
the sum of the rotational and translational kinetic energies: 



K=]^{mM) + '^{{R + r)H^ + z^ 



(6.16) 



Here, il denotes the angular velocity written in the body frame and I is the inertia tensor that is a 
diagonal 3x3 positive definite matrix I = diag{Ii, I2, I3) by our choice of body frame. 

Symmetries 

Since both the orientation and the height of the space frame is arbitrary, we expect to have a, G = 
M.X S^ symmetry. The group action on a point {g, z, {}) in the configuration space Q = S0{3) xM.x S^ 
is defined as: 

{z',^'):{g,z,^)^{m^,g,z + z',^ + ^') 

where 

/ COS??' -sint?' \ 
m^> = sin^' COS??' . (6.17) 

V 01/ 

This is a left action that with our choice of Euler angles is locally given by 

{z', ^') : (y^, 9, ij, z, ^) ^ {^ + ^' , 9,i:,z + z' , ^ + ^'). 

One can check that both the kinetic energy and the constraints are invariant under the lift of this 
action to TQ. 
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The action of G = M x S*^ on Q is free and proper and the shape space Q = Q/G = 5*0(3). In 
our local coordinates the orbit projection p : Q ^ SO (3) is given by 

p{<f,ij,e,z,^) = {^,^,e), (6.18) 

where ^ = ip — i), and {C,ip, 9) are Euler angles for 5*0(3) chosen with respect to the x-convention. 
The vertical subbundle Vp is spanned by 

{Odd 

Vp = span ijr + T^^ IT 
{ dip ov oz 

On the other hand, the 3- dimensional constraint distribution D is spanned by the G- invariant vector 
fields 

cos 6 suiiiD — 'd) d , „^ d sm.iip — 'd) d 



sm.9 dip 89 sin^ dip 

^d cos 9 cosiip — 'd) d . , ^. d cosiip — d) d 

y2 = R^ ^ -^ - si^^ - ^bw + n 7^' 6.19 

oz sin9 dip o9 smfc* dip 

d f R \ d 

^3 



dip \R + r J d'd 

It is then clear that Vp r\ D = {0} everywhere on Q. Therefore we have TQ = D ® Vp, the 
dimension assumption is satisfied and we are dealing with a Chaplygin system. The vector fields 
Yi,Y2,Y3 defined above verify the requirements of step 1 in the algorithm. Moreover, the condition 
in step 3 is vacuously satisfied. We thus concentrate in examining the condition in step 5. 

The projections Ya = Tp(Ya) are given by 

^ cos 9 sin C d d sin C d 

y\ = ^^r^T^ + cos^- ' 



sin 9 d^ d9 sin 9 dip ' 

-- cos 9 cos $, d . d cos C, d 

I2 = ^— — 7— - sm4— - + 



Y, 



3 



sin 9 dC, d9 sin 9 dip ' 

r \ d 



R + rJ d^ 
Following step 5, we consider the one- form Q that satisfies 

-^/<y^ 9 ( cos^sinA d , ,^,, d /sinA „ .^ „ ,x cos^cos^ , , ^- n ,\ 

^C^i = -^c ^1^ + 7m ^°^^ + 7^ -^ + h{i.9,ip) = — -^ + m,9,ip), 

04 \ smfc* / d9 dip \sin9J sm9 

^,:r>, d f cos^cosA d , . ,^,, d /cosA ^ ,^ ^ ,x cos^sinf ^ ,^ ^ ,x 

C^C^2 = wp —^ + 7^ -Sm + T7- — I +/2 e,^,^ = —^ + /2 e,^,^ , 

dE, \ sin9 J d9 dip \sin9 J sin9 

^{Ys) = ^ (]^) + /3(e,^,^) = M^,9,iP), 

where fa{^,9,ip) = G^^ for a = 1,2,3. The one-form u is uniquely determined by the above condi- 
tions (this is always the case when the dimension assumption holds). Using the explicit expressions 
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for Ya we find 

UJ = ( ) hd^+ i -—- + hcos^ - hsmn de + 

\ r J y smt^ / 

+ I sin 6^ sin ^/i + sin 6 cos ^/2 + I ) cos 9f^ J dif} . 

According to the algorithm, since (i^, 9, ip) are coordinates in a chart that is open and dense in 5*0(3), 
and the first de Rham cohomology group of 5*0(3) is trivial, a preserved measure exists if and only 
if u is exact. 

With the aid of a symbolic algebra software one can compute explicit expressions for fa(^,9,ip). 
These expressions are too long to be included in the present paper but we can provide the MAPLE"'"'^ 
file upon request. Below we discuss our findings. 

If the sphere is homogeneous (/i = I2 = I3) one gets /i = /2 = /s = which implies u = 
—dln{sm{9)) and there exists an invariant measure. This result is well known (see e.g. [B]). We will 
give an explicit formula for the preserved volume and the reduced equations of motion at the end of 
section 16.51 

In the general case, the condition that Q is exact implies that it is closed. This implies that the 
functions 

Gi(e,^,^):=(^)|^-|(/icose-/2sin0, 

G2(e,^,V^) := (^) ^ - ^ (sin^sine/i + sin^cos^/s + (^) cos^/3 

vanish identically in the chart < $ < 27r, 0<^<7r, 0<^< 27r. However we found that 
I TT n n\ f R \ I \/3mi?^ 

^1 ( o' o' T7 ) = (-'2 — h 



2'2'37 ' ' \R + rJ \A{{h + l2)mR^ + hl2 + {mR^f) 

Hence, a necessary condition for the existence of an invariant measure is that Ji = /2. Under this 
assumption one finds: 



^"2 77' T' o = (^3 - h 



.2' 4' 37 '" '' \{mR^{h + h) + 2hh){h + mR^) 

that can only vanish if Ji = I^. Therefore, we have shown 

Theorem 6.3. The reduced equations for a dynamically balanced sphere that rolls without slipping on 
the exterior of an infinite cylinder possess an invariant measure if and only if the ball is homogeneous 
(h = h = h). 

It is interesting to notice that in the limit when the radius r of the cylinder is infinitely larger 
than the radius of the sphere one recovers the classical Chaplygin sphere problem that possesses an 
invariant measure for arbitrary distributions of mass. The other limit case, when the radius of the 
cylinder r — > will be considered in the following section. 

To our knowledge, this is the first known mechanical example of a Chaplygin system that does not 
possess an invariant measure. A mathematical example had been given in [8]. 
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6.5. Dynamically balanced ball rolling on a vertical wire. Consider the motion of an inhonio- 
geneous sphere, whose center of mass coincides with its geometric center, that rolls without slipping 
on an infinitely long vertical wire. This corresponds to setting r = in the previous section. We 
follow the same notation and definitions for z and •&. 



Putting r 



0, the constraints of rolling without slipping (|6.15|) become 



^3, 



R{l>J\ sin -d — ijJ2 cos ^) . 



(6.20) 



Note that D is still generated by the invariant vector fields Yi, Y2, Y-^ defined in (6.19). However, the 



problem cannot be treated in the framework of the last section since the vector field Y3 becomes 

d 



Za- 



_ d_ 
dip d-d'' 



that lies on the vertical subbundle Vp. Therefore, in this case, the intersection V^p has constant 
rank 1 and we are no longer dealing with a Chaplygin system. Moreover, notice that 

dim(L) + Vp) = 4 < 5 = dim(TQ), 

at all points q E Q. Hence the dimension assumption is not satisfied. 

We follow the steps of the algorithm to analyze the existence of an invariant measure for the 
reduced equations. The first thing that we need to do, coming from step 1, is to redefine the vector 



fields Yi, Y2 in (6.19) so that they are S-orthogonal to Za=i- For the sake of clearness we change the 
notation of the vector fields in (|6.19) to 



Xa= 



X. 



cos 6 sin{(f — -t}) d 



dz 



sin 9 dip 

cos 9 cos{ip — -d) d 



„, d sin(o9 — {}) d 

cos y. -^)^+ '": n TTT. 
d9 sm 9 dip 



sin^ 



. , „, (9 cos(o9 — ■)!}) d 

sm(v9 - ?9)— H — -^—. 

dp d9 sin 9 dip 



and define: 



Yn '.— Xn 



S(Xa, Za=l] 



^a=l 1 



a 



1,2. 



S(^a=l,-^a=l) 

For simplicity we only consider the case in which the ball is axially symmetric. We assume that the 
inertia tensor is J = diag{Ii,Ii, J3). In this case the local expression for the metric tensor S is 

g = (Ji sin^ 9 + Is cos^ 9) dip"^ + h d9'^ + Jg dip"^ + Jg cos 9 dp> ® dip + mR^ d^"^ + mdz'^, 

and we get 

sin((y9 — {}) d 



Y 



a=l 



d9 



Y^ 



— {Is + 'mR'^) cos 9 sm{p — {}) \ d d 

sin 9{h sin^ 9 + I3 cos^ 9 + mR"^) J d^ ^ ^^^^"^ ~ '~ 
(Ji — /a) cos6'sin6'sin((y9 — 'i9)\ d 

d 



+ 



sin 9 dip 



+ 



4+ 

dz 



+ 



hsin'^ 9 + hcos^ 9 + mR^ 

— {I3 + niR^) cos6'cos(v9 — t9) 
sin 9{Ii sin^ 9 + I3 cos^ 9 + mR"^) ) dp} 
(/i — /s) cos6'sin6'cos(v9 — ?9)\ d 
Ji sin^ ^ + /a cos2 e + mi?2 ) M 



„, d cos(o9 — {f) d 

sm^ - ^ 7^ + ^^TT^TTT + 

d9 sm 9 dip 
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It is readily seen that the commutators 

[^a=i,l^J=0, a = 1,2. 

Therefore, the coefficients Cf^^^ « = ^^ ^^ (^i ^i ^^^d the condition in step 3 of the algorithm 
is verified. We now concentrate in the the study of the condition given in step 5. Recall that 



the principal bundle projection p : Q — )■ 5*0(3) is given by (6.18). The projected vector fields 
Ya = Tp{Ya) are then 

-- cos 6* sin ^ d d sin^ d -- cos 6* cos .^ d d cos^ d 

Ya=l = ^-^ 7^7 + COS^T-r + ^"^Try, ^0=2 = ^"^ 7T7 " Sm^TTTT + ^ 



sin 9 d^ 89 sin 9 dip ~ sin 9 d^ 89 sin 9 dip 

Following step 5 of the algorithm we consider the family 9^ of one-forms u that satisfy 
^,4^ , cos ^ cos f ^ /^ ^x ^,c? N cos 6^ sin f ^ ,^ ^^ 

a; y«=i = —?r^ + /i e,^ , wC^a=2 = — —^ + /2 e,^ , 

sm 9 sm 9 

where 

The independence of /i, /2 on ip is obvious since all objects involved in the calculation are independent 



of ip (see Remark 6.4 below). The above relationships imply that the family 11 can be described as: 

Q = a(e, 9) d9 + 6(^, 9) dip + X{d^ + cos 9 dip) , (6.21) 

where A = X{^,9,ip) is arbitrary, and 

cos 9 
a{^,0) ■= — ^^ + /i cos ^ - /2 sin ^, b{^,9) := /i sin 6^ sin ^ + /a sin 6* cos ^ . 
smt^ 

Step 5 of the algorithm states that a necessary condition for the existence of an invariant measure is 
that there exists a member of this family that is exact. 

With the aid of a symbolic algebra software (we can provide the MAPLE^^ file upon request), 
one gets: 

M^,9) = ^^^^^^^^^^cos^sin^cos$, M<l>,9) = "^"^'^'^ + "^f^^"' ' ''K os9sin9 sin<!> , 
where D is the determinant of the symmetric positive definite matrix 

9iZa=uZa=l) 

S(F„=i,i;=i) 9{Y^=i,Y^=2) 

S{Y^=2,Y^=i) S{Y^=2,Y^=2) 

The homogeneous case 

In this case we have Ji = I3 and hence /i = /a = 0. By taking A = we get 

cos 9 
Q = -^d9 =-dan(sm9)) 

8110.9 

which is exact. We shall prove that in this case there exists an invariant measure for the problem 
below by explicitly computing the reduced equations of motion and giving a formula for the measure. 
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The inhomogeneous case 



We now show that if Ji 7^ I3 then the one-form (6.21) is not closed for any smooth function A. 
This implies that no member of 91 is exact and, by step 5 of the algorithm, that there is no preserved 
measure. The condition du) = yields the relations 

dX da ,1 fdb dX\ db dX dX 

Combining the first two relations we get an explicit formula for A that is independent of ip. Substi- 
tution into the third relation gives the following necessary condition for the existence of an invariant 
measure 



db cos 9 / d'^b q^^^\ 



+ — ^1 T^TT^ + COS^,^) =0. 

However, using MAPLE"""^ we find 

db ^[^ J^\ =_^{ {h-h)\mB?f \ 

ae^sin^U^^^ °' ^e; ($=.,,=./4) 4 \h{I;^mR^){{mR^Y^hh^mR\h^h))' 

that can only vanish if I\ = I^. Therefore, we have shown that there is no preserved measure if 
h + h. 

Remark 6.4. The assumption that the ball is axially symmetric simplifies the analysis significantly 
since all the quantities that appear are independent of the angle i\}. This is a consequence of an 
additional symmetry of the problem that corresponds to rotations of the sphere about its axis of 
symmetry. The abelian Lie group G = Mx S^ x 5^ acts on the configuration space Q = 5*0(3) xMx 5*^ 
by the rule 



where the matrices 9^^/,9i^' are defined as in (6.17). This action is free and proper, and leaves 



the constraints and the kinetic energy invariant. Nevertheless, the methods described in this paper 
and in [42j cannot be used to analyze the existence of an invariant measure of the resulting reduced 
equations. The reason for this is that the rank of the intersection V^p = Vp H D is not constant 
throughout Q. If the axis of symmetry of the sphere is perpendicular to the vertical wire, then the 
rank of V^p is 2. In any other configuration this rank is 1. Therefore, we are forced to work with 
the action of the smaller symmetry group G = M. x S^ (for which the dimension assumption is not 
satisfied). o 

The reduced equations and the expression for the invariant measure for a homogeneous 
sphere that rolls without slipping on a vertical cylinder and wire 

The expression for the invariant measure for a homogeneous ball that rolls without slipping on a 
vertical cylinder or wire can be readily obtained as we now show. We treat both problems simulta- 
neously by considering r > where, as before, r denotes the radius of the cylinder (the case of the 
wire corresponds to the case r = 0). We begin by introducing the modified Poisson vectors 

OL := cos^g~^ei + sm.'dg~^e2-, (3 := —sm.'{}g~^ei + cos^g~^e2, j := g^^e^. 
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They form an orthonormal basis of M^ and can be considered as the columns of an element in 5*0(3). 
The constraints (6.15) can be written in terms of these vectors as 



7? 



R 



R + r 



(^,7) 



-R{n,f3). 



(6.22) 



Using the constraints one can obtain the following kinematical evolution equations for the Poisson 
vectors 



ex 



/3 



R 



R + r 
R 



R + r 

7 = 7 X f2. 



{j,n)i3 + (xxn, 
{-r,n)a + i3xn, 



(6.23) 



Recall that the Lagrangian of the system (6.16) is given by 

K=^-{m,n) + '^(^{R + r)^' + z' 

and, in view of the constraints ( 6.22[ ), we obtain the following equations of motion 

R 



in 



Iflxfl- \i 



R + r 



7 + A2i?/3 



(6.24) 



m{R + rYi^ = Ai, 
mz = A2, 



where Ai,A2 € M are Lagrange multipliers. By differentiating the constraints (6.22), we obtain the 
following expressions for the Lagrange multipliers: 

Ai = m{R + rf^ = mR{R + r){j,n), 
niR^ 



A, 



mz 



R + r 



{-r,n){a,n)-mR{(3,n). 



Substitution into (6.24) yields 



+ mR^E)n = inxn + mR^{a, n)a + 



niR^ 
R + r 



{j,n){ct,n)f3, 



(6.25) 



where E denotes the 3x3 identity matrix and we have used the orthonormality of ex., f3, 7. We can 
solve for fl in the above equation by obtaining an expression for (ck, f2) in terms of f2 and ey.,f3. 
Taking inner product on both sides of (6.25) with (J + mR^E)~^ex and isolating (ck, Q) gives 



(a, 17) 



iri X fl + fl^ (7, n) {ex, n)(3 , (l + mR^Ey^ex 
1 - mR^{{I + mR^E)-^ex, ex) 



Once the above expression is substituted into (6.25) and the system is complemented with equations 
(6.23), we obtain a closed system for f2,Q;,/3,7 that should be thought as variables in the reduced 
space that is obtained by eliminating d and z by the symmetry. 
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If the ball is homogeneous with inertia matrix I = I ■ E then Ji7 x 17 = and also {a, ft) = by 
orthonormality of ex and /3. The equations simplify to 

[K + rj[I + mn^j 
a = {j,n)f3 + axn, (6 26) 

7 = 7 X 17. 
A direct calculation that uses (7, /3) = (7, a) = 0, shows that, for any r > 0, the measure da Ad/3 A 



d^ A dQ, is preserved by the flow of (6.26). Therefore we have 



Theorem 6.5. The reduced equations for a dynamically balanced, axially symmetric sphere that rolls 
without slipping on an infinite vertical wire possess an invariant measure if and only if the hall is 
homogeneous. 

7. Conclusions and future work 

We have presented a geometric setup that allows us to obtain necessary and sufficient conditions 
for the existence of an invariant measure for certain types of nonholonomic mechanical systems with 
symmetry. Our methods have been successfully applied to prove the existence and non-existence of 
an invariant measure for concrete problems. 

Moreover, our geometric framework provides a setup that might be useful to determine conditions 
that guarantee the existence of an invariant measure for systems with a particular kind of geometric 
structure (such as LR systems [HIES])- A particular class of systems that we plan to study concerns 
nonholonomic systems for which the reduced configuration manifold is a homogeneous space. Some 
interesting results concerning the existence of a preserved measure for these systems have been given 

in [HEg. 

We believe that, following the research line of [21 [12], the methods developed in this paper can 
be extended to consider nonholonomic systems subject to affine constraints in the velocities. To 
our knowledge, the problem of the existence of invariant measures for such systems has not been 
considered so far. 

On the other hand, it would be interesting to exploit the generalized nonholonomic connection. 



introduced in Section 4.1, in order to study the geometry of symmetric nonholonomic systems which 



don't satisfy, in general, the dimension assumption. 

Appendix. Volume forms on vector bundles 



Let T : E -^ Q he an orientable vector bundle, over an orientable manifold Q. If a is a section of 
T : E -^ Q, we can define its vertical lift a"^, which is the vector field on E given by 

d_ 

dt\t=o^ 

If {e^} is a local basis of sections of E and a = a^e^ then 



"''(7g) = :tt (7g + ta{q)), for 7, G Eg. 



d 
dpi3 
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where p^ are the coordinates on the fibers of E obtained using the basis {e^}. 

Lemma A.l. [23] Let u be a volume form on Q and Q be a volume form on the fibers of E* . Then, 
there exists a unique volume form u Afl on E such that 

V A fi(Zi, . . . , Z^, a^, . . . , a^) = z/(Zi, . . . , Z„)fi(ai, . . . , a„), (A. 2) 

for ai, . . . , a„ G r(r) and Zi, . . . , Z^ vector fields on E which are r-projectable on the vector fields 
Zi,...,Zm on Q. 

Locally, if (g*) are local coordinates on an open subset U ^ Q and {ca} is a basis of sections of E* 
such that 

i^ = dq^ A--- A dq^, Q = d A ■ • ■ A e„, 
then 

u An = dq^ A--- Adq"" A dpiA--- A dp^. (A.3) 

A volume form $ on E' is said to be of basic type if 

£^v$ = 0, Va G r(r). (A.4) 



Using (A.3), it is easy to prove that the volume form i/ A f2 is of basic type. In fact, we have the 
following result 

Proposition A. 2. A volume form ^ on E is of basic type if and only if there exists a volume form 
u on Q and a volume form Q on the fibers of E* such that 

^ = u AQ. 

Proof. Suppose that $ is a volume form on E of basic type. 

Let z/q be an arbitrary volume form on Q and Qq a volume form on the fibers of E*. Then we can 
assume, without the loss of generality, that 

$ = e/z/o A fio, with a G C°°{E). 



Now, using (A.4), it follows that 

rfa(7") = 0, V7Gr(r), 
which implies that 5" is a basic function with respect to the vector bundle projection r : E ^ Q. In 
other words, there exists a G C°°(Q) such that a = a o t. 

Thus, if we take 

u = e"vQ, VL = Qq, 

we have that $ = z/ A fi. D 

Next, suppose that D is an orientable vector bundle over Q with rank n and vector bundle projec- 
tion T : D ^ Q and that S is a bundle metric on D. We define the section A^ of the vector bundle 
A^r : A"D — )■ Q that is characterized by the condition 

AS(g)(«J,...,<) = l, 

if {u^, . . . ,«"} is an orthonormal basis of D* with positive orientation and q E Q. Locally, for an 
arbitrary basis {cq,} of r(r), we have 

A^ = ei A • • ■ A e„. 

A/detg(ea,e/3) 
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The dual version of A^ is the section Q^s of maximal degree of D* that is characterized by the 
condition 

^As(g)K,---,0 = 1, 

if {ul, . . . ,u1} is an orthonormal basis of Dg with positive orientation and q E Q. It has the local 
expression 

^As = \/dd~g(^J|^^e^ A--- Ae**, (A.5) 

where {cq} is a local basis of r(r) and {e"} is the dual basis of r(r*). 

Next, we will apply the previous construction to the particular case when the vector bundle 
Td : D ^ Q is just the tangent bundle to Q. 

Let g he a Riemannian metric on an m dimensional orientable manifold Q (that is, g is a bundle 
metric on TQ). Then, we can consider the volume form Q/^g on Q. It is just the Riemannian volume 
form which we will denote by Ug. The local expression of Ug is 

Ug = Jdetg{Xi,Xj)a^ A ■ ■ ■ A a™, 

with {Xi} a local basis of vector fields in Q and {«*} the dual basis of 1-forms. 

Finally, suppose that t : D —^ Q is an orientable vector bundle endowed with a bundle metric S 
over an orientable Riemannian manifold {Q,g). Then, it is clear that 

$ = z/g A A^ 

is a basic volume form on the dual bundle r* : D* ^ Q to D. 
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